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Abstract 

We consider the discrete Boussinesq integrable system and the compatible set of differential 
difference, and partial differential equations. The latter not only encode the complete hierarchy 
of the Boussisesq equation, but also incorporate the hyperbolic Ernst equations for an Einstein- 
Maxwell- Weyl field in general relativity. We demonstrate a specific symmetry reduction of 
the partial differential equations, to a six-parameter, second order coupled system of ordinary 
differential equations, which is conjectured to be of Garnier type. 

1 Introduction 

One of the outstanding universal features of integrable nonlinear evolution equations is their ap- 
pearance together with a structured infinite hierarchy of commuting flows, i.e. an infinite set of 
mutually compatible partial differential equations (PDE). This feature reflects the extremely high 
degree of symmetry hidden in the system, and it has been exploited in many ways leading to the 
intimate connection between the associated hierarchies on the one hand, and the representation 
theory of infinite-dimensional Lie algebras and loop groups, on the other hand, cf PP, [2] and ref- 
erences therein. However, the point of view that one single equation in the associated hierarchy 
possesses a dominant role is misguided and, in fact, one should consider the infinite set of PDEs 
altogether, which constitutes the hierarchy, as being the integrable system. 

It is remarkable that there are situations in which the complete infinite hierarchy of flows of an 
integrable evolution equation, can be encoded in a single PDE. A prime example of such an equation 
was given recently in [5], and it is represented be a fourth order scalar PDE. The complete hierarchy 
of the Korteweg-de Vries (KdV) equation follows by systematic expansions from this PDE. In the 
present work we consider another hierarchy generating PDE, namely the one that encodes the 
complete hierarchy of the Boussinesq (BSQ) equation. This will be achieved by deriving, first, the 
compatible discrete integrable system. 

The construction of the latter is based on the infinite-matrix structures, introduced in [I]. In 
particular we consider those structures which serve as a representation of the loop algebra built from 
the hierarchy of flows, of the second member in the Gel'fand-Dikii hierarchy, namely the Boussinesq 
equation [5J • The basic ingredient consists of a linear integral equation for (infinite) basic functions 
depending on a complex parameter and on the coordinates of the system. The coordinates can be 
chosen to be discrete as well as continuous, and exactly this freedom in the coordinate choice, is 
responsible for making the integral equations such a convenient tool to derive discrete integrable 
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systems. This is the subject matter of section 2, which ends up with a parameter family of discrete 
equations, namely the lattice BSQ equation. 

In section 3, the multi-dimensional consistency of the lattice BSQ equation is tested, by em- 
bedding the two dimensional lattice into a higher dimensional lattice. This aspect of the discussion 
forms a completion to the study of the lattice BSQ integrability as performed in [£]. The fact 
that the evolution in the extended three dimensional lattice is well posed reveals, in an essentially 
algorithmic fashion, the Lax pair of the lattice BSQ equation. Section 4, proceeds with the deriva- 
tion of the compatible set of differential-difference equations of the lattice BSQ equation. This will 
be achieved by imposing on the dependent variables an additional dependence on the continuous 
lattice parameters, in such a way that certain compatibility conditions are satisfied. 

Sections 5 and 6, deal with the continuous compatible PDEs and ODEs, respectively. First, we 
present the generating PDE of the Boussinesq hierarchy, consisting of a two-parameter family of 
coupled fourth order system of PDEs, in two dependent and two independent variables. Next, we 
illustrate the connection of the latter PDEs with the Boussinesq hierarchy, using the Lax pair of 
the relevant PDEs, which was derived in the preceding sections. A prominent position, among the 
subsystems which are included into the solution space of fourth order parameter family of PDEs, 
is possessed by a second order system. It is the heart of the Einstein's field equations for plane 
symmetric spacetimes, in the presence of certain massless fields. In the context of general relativity, 
the latter equations are known as the Ernst equations for an Einstein-Maxwell- Weyl field. Finally, 
we consider a particular type of similarity solutions of the system of PDEs under consideration, 
which are built from the solutions of a second order coupled system of ODEs, involving six free 
parameters. The latter may be considered to be associated with the hierarchies of the Painleve VI 
equations, 

The paper concludes with perspectives, where various avenues for future research are discussed. 

2 The discrete Boussinesq system 

In order to derive the lattice Boussinesq (BSQ) system, we shall use a linear integral equation of 
the following form 



cf [3]. The vector u(k), has i € Z entries, each of them taking values in the complex numbers, and r 
is a vector with entries given by n = (i 6 Z). The variables £ and k are complex, and u is a iVth 
root of unity, i.e. uj = exp(2i-7r /N). The integrations in Q are performed over an arbitrary contour 
C in the complex A:-plane, and over an arbitrary measure dX(£). In principle, the only restrictions 
on the measure and the contour consists of the requirement of the existence and uniqueness of the 
solution of the integral equation (^Q), once the inhomogeneous term is specified. Along with the 
infinite vector u, which solves the integral equation lfl|). we can introduce an infinite-component 
matrix potential U, given by 



where r denotes the transposed vector of r. The aim is to derive for the matrix potential U, 
as well as for the vector potential u, a set of (infinite-) matrix equations, by imposing on these 
objects a dependence on additional variables and certain transformation rules. The variables can 
be continuous as well as discrete. Adopting the latter case, we consider two operations, namely the 
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maps 



u(n, m, k) i — > u(n, m, k) = u(n + 1, m, k), u(n, m, k) i— > u(n, m, k) = u(n, m + 1, k) . (3) 

These maps have the interpretation of shifts in a two-dimensional lattice, but equally well of 
Backlund transformations generated by a transformation of the complex function p(k), which is 
not specified yet. The transformation rules which are imposed on p{k) are 

p(k; p, q) i-> p(k; p,q)= P + - p(k; p, q) , p(k; p, q) i-> p(k; p,q) = g - p(k; p, q) , (4) 

p + ojk q + ojk 

where p and q are the complex lattice parameters, associated with the lattice variables n and m, 
respectively. Using the linear integral equation we may determine how u changes, under the 
action of the transformations (@J. Restricting to the special case N = 3, i.e. lo = exp(2i-7r/3), the 
equations which determine the change of u, are given by 

(p + uk)u = (p + A-UO)u, (5a) 
(p + k){p + u) 2 k)u = (p + ojA) {p + J 2 A) u 

-U (luO{p + lu 2 A) +lu 2 (p + lu 2 *A)0) u, (5b) 

together with the same equations by replacing ~ and p, with ~ and q, respectively. The matrices 
A and *A, define the operations of index-raising when multiplied from the left and from the right, 
respectively, i.e. (Ait)j = iti+j. The matrix O is a projection matrix which singles out the central 
element of the object which upon acts, i.e. (Ou)i = uo5io, where <5y the Kronecker delta symbol. 
Integrating (JSJ) over the contour C and the same measure dA(£), we obtain the following nonlinear 
system of equations for the matrix potential U 

U{p + u *A) = (p + A — UO)U, (6a) 
Uip + ^ip + u 2 f A) = (p + cuA)(p + uj 2 A)U 

-U (uO(p + uj 2 A) +lu 2 (p + u 2 *A)0) C7. (6b) 

Using the identity a; 2 + c<j + 1 = 0, the above equations may be rewritten as follows 

Uip + u^) = (p + A-UO)U , (7a) 
U (p 2 - up f A + u? (*A) 2 ) = (p 2 -pA + A 2 )C/ + (7 (pO-OA-w *AO) U. (7b) 

From the above infinite set of relations for the components of the potential U, we restrict ourselves 
to the equations evolving its (0,0), (1,0) components from eq. (f7a|) and the (0,0) component from 
eq. (|7b|) . together with the same equations by interchanging ~ and p, with ~ and q, respectively. 
A simple elimination process, using the discrete operations, leads to a closed system of equations 
for the variables (Uqq, U\q, Uqi). The omitted latter set of discrete equations, obtains a simple form 



by introducing new dependent variables (u, v, w) by the relations 

u = Uqo - (np + mq) , (8a) 

v = U\q — (np + mq)Uoo + \n(n + l)p 2 + nmpq + ^m(m + l)q 2 , (8b) 

w = —luUqi — (np + mq)UoQ + ^n(n — l)p 2 + nmpq + ^m(m — l)q 2 . (8c) 
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In terms of these new variables, the aforementioned set of partial difference equations (PAE), reads 



w 
w 

w 



UU — V 
UU — V 



uu 



v + 



t — a 



u — u 



(9a) 
(9b) 

(9c) 



where r = p 3 , a 
system. 



g 3 . The system of PAEs ©, will be referred in the following as the lattice BSQ 



3 Three dimensional consistency of the lattice Boussinesq 

The lattice BSQ system (|9*|). represents a compatible parameter family of PAEs, in the sense that it 
can be embedded consistently in a multidimensional lattice, on which the evolution is well posed, 
0. Recently, this property has been successfully exploited for quadrilateral scalar lattice equations 
to obtain a classification of the latter, [S]. 

In order to illustrate the above statement, we extend the lattice into a third dimension by 
introducing a new lattice variable I associated with a new shift, which will be denoted by T , and 
a new complex parameter Embedding the system of equations © into the generated three 
dimensional lattice, we get 



w 



w 



uu — v, 



uu — v, 



uu — v, 
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uu — v + 
uu — v + 
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(10a) 
(10b) 
(10c) 



Using from (|10j) . those equations which involve the shift of the variable w in one direction, and 
taking into account the relevant compatibility conditions, we can determine how u and w, are 
shifted in two directions. The corresponding equations for the variable u read 



v — v 



u — u 



— v — V 
u = = — — 

u — u 



— v — V 

u = — — — 

u — u 



(11) 




(u,v) & 



(u,v,w) 



Figure 1: Three dimensional consistency of the lattice Boussinesq system. 
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Let us now impose the following initial values 



(u,v,w), (u,v), (u,v), (u,v), (12) 

at the vertices depicted by • in Figure ^ Successively, using equations lfTU|) . (|11() we may evaluate 
the updates 

(u,v), (u,v), (u,v), (13) 

which are assigned to the vertices depicted with o. Then, from these data and the evolution of the 
lattice equations (|1()|) . (|11[) we are free to calculate the values at the final point, in three different 
ways. In principle, there is no particular reason for these three values to coincide, unless the 
PAE's © enjoy this special property. The confirmation of this test is crucial in that it verifies, at 
first stage, the consistency of the embedding of the lattice BSQ into a multidimensional lattice. By 
simple and straightforward algebraic calculations, and eliminating all terms shifted in two directions 
in favour of terms shifted in one direction, one finds that 

S = M+ K-ffi-^-K-;><g-;> (14 a) 
(v — v)[u — u) — [v — v){u — u) 

i = r>+ ( ,l"^g"5 g "g"^"5 g , (14b) 
[v — v)(u — u) — [v — v)(u — u) 

I = M+ K-ffi-^-K-;»<g-;; , (14c) 

(v — v){u — u) — (v — v)(u — u) 

independently of the way of calculating the values at this point. This fact can be easily verified, 
since the fractional terms in equations (|14|) . remain invariant under any permutation of the three 
lattice shifts, along with the corresponding lattice parameters. 

Having established the three dimensional consistency of the lattice BSQ (using no more infor- 
mation than the equations themselves), it is now straightforward to derive the linear representation 
for the lattice BSQ system. This is accomplished by appointing (u, v) as projective variables, and 
identifying them with the ratio of homogeneous variables (ipo^i,^), as follows 

u=—, v = — . (15) 
Wo Wo 

Substituting (u, v) given by equations (|15|) , into the relevant equations fi()[) , (|11|) and setting 

fa = -UljjQ lj)Q = -uip Q + Vi , (16) 

we obtain 

ipi = -vy -ip2, $2 = {t — C + w u — uv)^o — wipx + U1p2 , (17a) 

ipi = -vipQ + W2, W2 = (f - C + wu - uv)ipo - wWi + UW2 ■ (17b) 

The linear system (|16[). (|17p can be written in a matrix form as follows 

-u 

x/> = Li/)= \ -v 1 | i/>, ip = MW = ( -v 1 | V , (18) 



—u 


1 


°\ 




—v 





1 W, 
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— w 


u ) 
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where tp = *(^cb V'l; V^)) and the entries of the matrices L and M which are denoted by > and <, 
are uniquely determined by the relations 

det L = t - (, det M = a — ( , (19) 

respectively. It should be noticed that the compatibility condition, LM = ML, of the linear system 
(|18|) for every value of £, holds on weaker equations than the lattice BSQ system 0. The validity 
of the previous compatibility condition on the solutions of the lattice BSQ, implies that the linear 
difference system (|18|). can be solved for every fixed value of £. Thus, we may replace the column 
vector ijj, by a matrix with values in GL(3, C), and which is formed by putting together three 
independent solutions of the system (j!8j) . 



4 The differential-difference compatible system 

We now impose on and subsequently on the dependent variables, an additional dependence on 
the continuous lattice parameters (r, a). This will be done by imposing on the linear system pHj). 
certain first order linear differential relations for *S>, in such a way that the following compatibility 
conditions 

= d T ^ , 9 CT § = d^> , (20) 

are satisfied. Partial differentiation will be denoted by a comma, or by <9, followed by the variable(s) 
with respect to which the differentiation has been performed. From the linear system (|18|) . where 
now if) is replaced by the matrix \E', it follows that 

det* = (r-C)>-C) m , (21) 

up to a non significant constant which may be fixed to unit, without losing generality. The previous 
observation and the matrix identity 

*(,,.-',= (22) 
suggest that, a potential simple form of the first order linear differential equations for ^, is 

*t = ^7 * , «V = -^7 * • (23) 

Here, ^4 and i? take values in gl(3, C), they do not depend on £ and 

tvA = n, trB = m. (24) 

The compatibility condition d T ^ = d^, for every value of £, applied on the relevant equations 
(|18|) . (|2cS|) . leads to the following system 

G :T + A + AA- AA = 0, (25a) 

GA-AG = 0, (25b) 

where 

G = L-(t-()A, (26) 
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and the matrix A is given by 

4=10 01. (27) 




After a lengthy elimination process applied to the system of equations (|25|). using the discrete 
operation, and taking also into account equation (|9a|) and the relevant trace relation from (|24|) . 
one obtains the explicit form of the matrix A, and simultaneously the compatible system of DAEs 
for the dependent variables. It turns out that the compatibility conditions (|25|) . amount exactly to 
the equations obtained by identifying two different dyadic forms for the matrix A. Specifically, A 
is given by 

Aij = cii a'j = Ci c'j , (28) 

where 

ai = (1, u, v) , a'i = (n + vu jT — uw^ T , w tT , —u jT ) , (29a) 

Cj = (— u T , — v )T , n + wu :T — uu T ) , c\ = (w, —u, 1) , (29b) 

i = 0,1,2, and j__ denotes backward shift in the lattice variable n. The identification of the two 
dyadic forms of A, leads to the following differential-difference equations (DAE) 

_ _ „ n 

» T = M« T , W jT = UU^ T , V + W — U U H =0. (30) 

The variables v and w can be easily eliminated from equations (|3U|) . leading to the following scalar 
DAE for the variable u 

— ) = (u — u) u + (u — u)u T . (31) 

V/.T ~ ' 

Since by construction the system is symmetric with respect to the lattice variables and the associ- 
ated parameters, the replacement f, n,r) with (^,m,a) in equations (|2U)) - H31|) . delivers the matrix 
B and the associated DAE for the variables (u,v,w). Explicitly, the matrix B takes the dyadic 
expressions 

B = h % = di d) , (32) 

where 

bi = (l,u,v), b'i = (m + vu )CT - uw y „ , w tCT , -u }(T ) , (33a) 

di = (-u o- , -v t(J , m + wu )CJ - uv )Cr ) , d'i = (w, -u, 1) , (33b) 

i = 0, 1,2, and^ denotes backward shift in the lattice variable m. The system of the compatible 
DAEs is 

nrn 

v CT = uu a , w CT = uu a , v + w — uu-\ = 0, (34) 

which by eliminating the variables v and w, leads to the following scalar DAE for u 



in 



u, 



(u — u) u a + (u — u) u iCT . (35) 



The DAEs in equations ()30|) and (|34[1. represent the Zietero-Backund transformations between 
copies of the compatible set of PDEs, which we derive in the next section, for different values of 
the parameters n and m. This issue, along with the derivation of auto-Backlund transformations 
of the continuous systems will be the subject of a future publication. 
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5 Compatible partial differential equations 



We now focus on the dependence of ^f, and the dependent variables, on the continuous lattice 
parameters (r, a) S C 2 . Since we drop for the moment the discrete operations, the variables which 
are assigned to different lattice vertices are interpreted as new dependent variables. However, 
we point out that, we may retain the different lattice labels on the associated fields, since by 
construction the whole set the PAE's, DAE's and PDE's are mutually compatible to each other. 



5.1 The generating PDE of the Boussinesq hierarchy 

The compatibility condition ^ jTa = ^,ar on the linear system (|23j) . for every value of £, leads to 
the following system of matrix PDEs 

(a - r)A ff + [A , B) = , (a- t)B :T + [A , B] = . (36) 



Let us now choose the dyadic forms of A and B, given by (|29a|) and (|33a|) respectively, and which 
for convenience we rewrite below, along with the relevant trace relations given by equations (|24|) . 

Aij = aid'j , a 1 a'i = n , (37a) 

B tj = hb'j , 6* b\ = m , (37b) 

i,j = 0,1,2. Indices are raised and lowered using the Kronecker symbol 5{j, and summation over 
repeated upper and lower indices is understood. From equations (|36|) . it follows that we may 
consider a more general case where n and m are parameter functions of r and a, respectively. In 
the following, we adopt the latter option, unless otherwise stated. We shall only require that when 
the relevant discrete operations are performed, n and m, vary by unit steps only. 

Inserting A and B, given by Q37JI into equations (|36|) . and after simple algebraic manipulations, 
taking into account the fact that ao = &o = 1> we arrive at the following system of PDE's 

(r — a) a^c — a? bj (a, — b{) = 0, (r — a) a' i a + a? b'j a[ — Va'j b\ = , (38a) 

(r - a) b i>T - b> a'j (oj - h) = 0, (r - a)b' ir + a j b'j a- - b> a'jb'i = , (38b) 

i,j = 0, 1, 2. From the right set of the above of equations we find that a' ia = b[ T . Using this fact, 
and the explicit form of a! i and b\ , given by (J29a|) and (J33a|) respectively, we introduce the variables 
(/o, fi, h)> by the relations 

fi=w, f 2 = -u, fo, T = n- a 1 f iiT , f 0>a = m-b % f iy(T , (39) 

i = 1,2, and which satisfy 

< = fi,r, b'i = f ha , (40) 
i = 0, 1, 2. In terms of the new variables fi, the system of equations (|38[) reads 

(r - a)a i>a - (m + (a j - b?) f j>ar ) (oj -64) = , (41a) 
(r - <7)6 i)T - (n - (a? - V) f hT ) (at ~h) = , (41b) 
(t - cr)fi, T a + mfi,r ~ nfi,o- + Cij(a j - b>) = 0, (41c) 

where i,j = 1,2. The components of the matrix (Cij), which is in general a rank two matrix, are 
defined as follows 

Cij fi,u fj,v fi,v fj,u j (^9 
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i,j = 1,2. Solving the linear system (|41c[) for the differences (a^ — 6j), and using the remaining 
field equations (|41a|) . ()41b|) . we obtain a fourth order coupled system of PDE's for the variables 
(/i, f 2 ), which is omitted because of its length. The latter system of PDEs are the Euler-Lagrange 
equations for the variational problem, associated with the Lagrangian density 

FE F E 

L = (t — a) h m(a) - + n(r) - . (43) 

Here, the scalar quantities F, E and J are given by 

F = f >T(T xf ;T , E = f >ra xf ta , J = f T xf >a , (44) 

where / is the vector with components (fi,f 2 ). In the following, we denote the parameter family 
of fourth order PDEs for the variables (/1, 72), by E(L) = 0. We conclude this section by noticing 
that the Lagrangian L remains invariant, modulo L and null Lagrangians, under the action of the 
projective group of linear fractional transformations on (/1, 72), i.e. 

(f f) ( Anfi+A 12 f2 + Ai 3 A 21 h + A 22 f 2 + A 23 \ , . 

UhM » y Asifi + Amh + ^ 33 > Asifi + A32h + Ass J » I ) 

where (Aij) € SL(3, C). If in addition the parameter functions n and m are constants, then L 
remains invariant under the affine base transformations 

(r, a) (e\T + e 2 , E\cj + e 2 ) , ei, e 2 G C, E\ / . (46) 

The symmetry transformations of L, given by (|45|) and (|46|) . fall into the class of the so-called 
divergence symmetries of the associated variational problem. Every symmetry of this type is 
inherited as a Lie-point symmetry, to the corresponding Euler-Lagrange equations, see |1U| for a 
detailed account on these topics. 

5.2 The Boussinesq hierarchy 

The coupled system of PDEs E(L) = 0, encodes the complete hierarchy of the Boussinesq equations. 
This can be illustrated, by imposing on ^ and on the associated dependent variables, an additional 
dependence on an infinite set of indeterminates x l = (x 1 ,^ 2 , . . .). Following the construction given 
in |12j . the dependence is given explicitly by certain (infinite) differential relations, in such a 
way that all these relations are mutually consistent to each other. To this end, we introduce the 
matrix H with components 

H-ij,T 0S( f j t r 1 Hij,cr fj,cr > (^^) 

i, j = 0, 1, 2, and the matrix Z by the relation 

1 \ 

(48) 

The existence of the matrix H for the variables (r, a), may be easily verified by using the system 
of PDE's (|41|) . In terms of the above matrices, the linear differential relations for which can be 
simultaneously imposed on the linear system (|23[). are given by 

8^ = (Z + [ d c Z,H])V, (49a) 
d 2 ^ = (Z 2 + [d^Z^H])^, (49b) 
8 i+3 * = -({di* + d % H)^ : (49c) 
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where 8i = 8 x i and i G Z/3Z. 

Taking the compatibility condition 8182^ = <9 2 <9i v I / , on equations (|49a|) and (|49b|) . we obtain 
a set of PDE's for the various entries of the matrix H. A straightforward differentiation and 
elimination process applied to the latter system, and identifying the main field variables (/i,/2) 
using (|47|) . we end up with the following set of PDE's 

cWi = 2(5 2 ) 2 / 2 -3a 1 (a 1 / 2 ) 2 -i(9 1 ) 2 a 2 / 2 + i5 1 a 2 (/|) + i(9 1 ) 4 /2, (50a) 
d x d x h = \8 x d 2 h + ^ 1 (a 1 / 2 ) 2 -i(5 1 ) 3 / 2 . (50b) 

Then, the final compatibility condition (<9i) 2 $ 2 / = <9 2 (<9i) 2 / on the above system, leads to the 
Boussinesq equation 

3(5 2 ) 2 Q + (8 x fQ + 6(d!) 2 (Q 2 ) = , (51) 

where Q = —81/2- Similarly, higher members of the Boussinesq hierarchy can be derived from the 
system (|49j) . in a systematic way. 

5.3 The Ernst equations for an Einstein-Maxwell- Weyl field 

Now we discuss various subsystems which are included into the equations E(L) = 0. Such sub- 
systems arise, by imposing additional algebraic or differential constraints on the space of the in- 
dependent and the dependent variables and on the parameter space, in a consistent way with the 
equations E(L) = 0, or equivalently, with the auxiliary system (|4*T|) . 

Assuming first that the independent variables are real, one such consistent algebraic constraint 
is given by the following set of relations on the auxiliary variables, 

((J 7~) y* (J T 

ai = h+ ° / 1 , a 2 = b 2 + - J 7 I , (52) 

J2 + h + / 1 / 1 h + J 2 + h / 1 

and on the parameters 

n + n* + 1 = 0, m + m* + l = 0, (53) 

where * denotes complex conjugation. Imposing the constraints (|52j) and H53|) on the auxiliary 
system (|4"T|) , the latter reduces to a second order coupled system of PDEs for the main dependent 
variables (/i,/ 2 ) solely, which is given by equations (|41c|) . Relabeling the variables (/i,/ 2 ) by 
(<P, £), and using the language of forms, the latter set of PDEs can be written as 

(£ + S* + $$*) (&{q * d£) - i A d£) = 2g(d£ + $*d$) A *d£, (54a) 
(£ + £* + (d(g * d#) - i /3 A d#) = 2^(d^ + $*d$) A *d<^ . (54b) 

Here, * denotes the Hodge duality operator acting on the basis of one forms by 

*dr = dr, *da = —da , (55) 

while the function g, and the one-form f3 are defined as 

Q=-(a-T), (3 = u(r)dT + fi(a)da , (56) 

where v{t) and n(cr), are real parameter functions of the indicated arguments. The system of 
PDEs (|54|) forms the basis of the Einstein's field equations for plane symmetric spacetimes, in the 
presence of a source-free Maxwell field and a Weyl neutrino field, and they are known in general 
relativity as the Ernst equations for an Einstein-Maxwell- Weyl field, |13j-|15j. 
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6 Symmetry reduction to a second order system of ODEs 



We are now interested in finding particular similarity solutions of the equations E(L) = 0. This 
will be achieved by, first, fixing the parameters functions n, m to be constants, and then by 
choosing a specific one-dimensional subgroup from the full group of Lie-point symmetries (|45[). 
(J46|). The requirement of invariance of the solutions under the one-dimensional subgroup, leads to 
the corresponding similarity solutions. 

Our motivation in finding such solutions is two-fold. First, one hopes to discover some explicit 
solutions of the fourth order coupled system E(L) = 0, which would provide physically significant 
solutions, in connection with the Ernst equations. In the present work however, we are interested 
mostly in writing down the system of the reduced ordinary differential equations (ODE), since the 
latter potentially define higher order analogous of the Painleve ODE's. 

To be more specific, we shall be interested in finding solutions of the equations E(L) = 0, which 
remain invariant along the orbits generated by the following vector field 1 

X = rd T + ad a + 2 ai fid fl + 2a 2 f 2 d f2 , (57) 

where a\, a 2 , arbitrary complex parameters. Requiring invariance, means that the function (f\, f 2 ) 
should satisfy the equations E(L) = 0, and also the compatible differential constraints 

Tfi,T + crfi,v = 2a i f i t = 1,2, (58) 

The general solution of the latter PDEs is 

T 

fi(r, a) = Fi (t) (ra) ai , t = — , i = 1, 2 . (no summation) (59) 

a 

Inserting the main variables in the form given by equation (|59j) . into equations E(L) = 0, we 
obtain a fourth order system of ODEs for the functions Fi, which is omitted because of its length. 
However, an effective way in dealing with the latter system of ODEs, is to consider the auxiliary 
system (j41j) . To this end, one needs first to find the vertical components of the vector field X, 
which are associated with the auxiliary variables aj, 6j, and subsequently the consistent invariant 
form of the variables a^, 6j. The corresponding symmetry generator X for the equations (|41j) . is 

X = X + (1 - 2a 1 )(a 1 d ai + hd bl ) + (1 - 2a 2 )(a 2 d a2 + b 2 8 b2 ) , (60) 

which leads us to conclude that the invariant form of the auxiliary potentials is 

ai (T,a) = P l (t)(Ta) 1 / 2 ~ a \ t . 10 , \ 

^J = Ri \J ) \Jy / 2^ t = ~, -1,2- (—.on) (61) 

Substituting the dependent variables a,, bi and fi given by (J6*T|) and (jSHJ), respectively, into equations 
(|41l) . we obtain the following system of ODEs 



-(t - 1) (tP! - (| - ai)P % ) = (m- tV 2 (P j - R j )(tF; - ajF,)) (Pi - Ri) , (62a) 

(t - 1) [tR'i + (i - cti)Ri) = (nt - t l ' 2 {P j - R j )(tF^ + ajFj^j (Pi - Ri) , (62b) 

t 2 (l - t)F" + t((n -l)t + m + l)Fl - ati((n - an)t - (m - (Xi))Fi = 

2t 1/2 (t 2 FlFj - aiFiajFj) (P j - R j ) , (62c) 



The symmetry generator H57II . generalizes the one used in |16|. for the symmetry reduction of a scalar 
subsystem, included into equations E(L) = 0, to the full Painleve VI equation. 
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where prime denotes differentiation with respect to t. Employing a general quadratic Ansatz on 
the first order derivatives of i^'s, we find that the system (|62[). admits two first integrals. Denoting 
the complex constants of integration by s± and s 2 , these first integrals are given by 



{Pi - Ri 



) 2 (t 2 F[ 2 - a 1 2 F 1 2 ) + (P 1 - /?!)(P 2 - R 2 ){t 2 F[F! 1 - a 1 a 2 F 1 F 2 ) 

F i + y Fi ) ( (2ai " " " " ^ ~ T Fl ^> ( (2ai " 1)Rl + n(Pl " 



+ - a 2 Fi'F 2 ) 



(Pxi? 2 - P 2 Pi) ai + a2 1 + (R1R2 - P1P2) 

Oi\ — OL 2 



Si , 



(63) 



together with the same relation, obtaind from (|63|) by interchanging the subscripts 1 «-> 2. Through 
a straightforward elimination and differentiation process applied to the system (|62[). taking into 
account the first integrals and introducing new dependent variables (0i,0 2 ) by the substitutions, 



&(t) 



an Fi(t) 



1,2. 



(64) 



we end up with a second order, coupled system of ODEs for the variables (Qi, 2 ), which involves 
six free parameters, namely (n,m, a±, a 2 , s\, s 2 ). In order to give an explicit form of the latter 
system of ODEs, we introduce the auxiliary quantities (Qi,Q 2 ) as follows 



Qi 



m + nt 



+ 



m-nt)g 2 (t- l)a 2 0i0 2 (l - 2 2 ) (t - \){a x - a 2 )(G 2 2 - 1) 



+ 



2a 1 (g 1 -g 2 ) 2a l {g 1 -g 2 ) 2a 1 {g l -g 2 ) 2 2a 1 (g 1 -g 2 ) 2 

(t - \)g x 2 {g 2 2 - 1) t(t-i){g 2 2 -i)g[ t(t- -i)g' 2 



+ 



2(Gi - g 2 f 



+ 



(65) 



and Q 2 is obtained from Qi by interchanging the subscripts 1 
system of ODEs is given by 



2a 1 (g 1 -g 2 ) 2 2a 1 (g 1 -g 2 y 

2. Then, the aforementioned 



atx a 2 
a\ — a 2 

t 

O 



(0i - g 2 )(QiQ 2 - Q[Q 2 ) + ((* - - (t + l))Qi 



(0i 



40i 



t-1 



Qi 



a.\ a 2 



(g 2 



402 ai(0i-0 2 ) s 



t-1 



) 2 - 40i + — 0i - 1) + — ^- Q 

ai ai (i — 1 



ai — a 2 



Q1Q2 



(66) 



together with the same equation obtained from equation (|66|) . by making the usual interchanges in 
the subscripts. 



7 Perspectives 

We finish our discussion by noting that, on the continuous level, the generating PDE of the Boussi- 
nesq hierarchy, namely the equations E(L) = 0, include the one of the KdV hierarchy. This is 
achieved by considering the case, where the dependent variables (/1, f 2 ) coincide, i.e. f ■= fi = f 2 - 
The generating PDE of the KdV hierarchy was given in [3] , and is the Euler-Lagrange equation for 
the variational problem associated with the Lagrangian 

r I \(f,ra) 1 / 2 f,r 2 f,(r\ fcrr\ 

C = {t-<j) + to— + n— . (67) 

J,rJ,a T-U \ } )Cr J jT J 
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and which we denote by E(C) = 0. The Lagrangian C can be reduced from the Lagrangian L 
given by equation (|43|). considering the aforementioned case, in a consistent way with the reduced 
equations. In [3], cf ^)], it was shown that a specific symmetry reduction of the equations E(C) = 
leads to invariant solutions, built from the full Painleve VI equation. Thus, we expect that there 
exists an intimate connection of the six-parameter system of coupled ODEs, presented in Section 
7, with the hierarchies of the Painleve VI and the Gamier systems |17| . The latter issue along with 
the isomonodromy problem of the relevant ODEs, will be left to a future study. 

The derivation of generating PDEs, which, in turn, provide a mechanism for obtaining new 
integrable systems, is a project of equal importance. To this end, the classification of scalar discrete 
integrable systems, obtained recently in [SJ, could serve as a starting point. However, some of the 
extra restrictions imposed on the classification scheme in would exclude, in a strict sense, the 
discrete BSQ system ©. Thus, the classification of quadrilateral integrable discrete systems, in 
more than one dependent variable, seems to be a more complicated problem. 

In the light of the intimate connection between the PDE's E(L) = and E(C) = 0, and the 
Ernst equations, the former fourth order integrable systems acquire a certain physical importance. 
This stems from the fact that the solution space of the hyperbolic Ernst equations is embedded into 
the solution space of the fourth order systems. Consequently, it would be interesting to consider 
to what extent, solutions of the larger continuous systems and their discrete compatible systems, 
could provide new exact solutions of the Einstein's equations, for various physically significant 
problems. Such a problem is the head-on collision of plane fronted gravitational waves, coupled 
with electromagnetic and neutrino waves, see ^S] and references therein. 

Acknowledgements 

The work of A.T. was supported by the European postdoctoral fellowship Marie Curie, contract 
No HPMF-CT-2002-01639. 

References 

[1] V. G. Kac, Infinite dimensional Lie algebras, (Cambridge University Press, Cambridge , 1990), 
third edition. 

[2] T. Miwa, M. Jimbo and E. Date, Solitons: Differential equations, symmetries and infinite 
dimensional algebras, (Cambridge University Press, Cambridge, 2000). 

[3] F. Nijhoff, A. Hone and N. Joshi, On a Schwarzian PDE associated with the KdV, hierarchy 
Phys. Lett. A 267 (2000), 147-156. 

[4] F.W. Nijhoff, G.R.W. Quispel, J. Van der Linden and H.W. Capel, On some linear integral 
equations generating solutions of nonlinear partial differential equations, Physica 119 A (1983), 
101-142. 

[5] F.W. Nijhoff, V.G. Papageorgiou, H.W. Capel and G.R.W. Quispel, The lattice GePfand-Dikii 
hierarchy, Inverse Problems 8 (1992), 597-651. 

[6] A.J. Walker, Similarity reductions and integrable lattice equations, (Leeds University, 2001), 
PhD thesis. 



13 



[7] F.W. Nijhoff and A.J. Walker, The discrete and continuous Painleve VI hierarchy and the 
Gamier systems, Glasgow Math. J. 43 A (2001), 109-123. 

[8] V.E. Adler, A.I. Bobenko and Yu.B. Suris, Classification of integrable equations on quad- 
graphs. The consistency approach, Commun. Math. Phys. 233 (2003), 513-543. 

[9] A. Tongas and F. Nijhoff, Generalized hyperbolic Ernst equations for an Einstein-Maxwell- 
Weyl field, |nlin.SI/0401010| preprint. 

[10] P.J. Olver, Equivalence, Invariants, and Symmetry, (Cambridge University Press, Cambridge, 
1995). 

[11] F.W. Nijhoff, Linear integral transformations and hierarchies of integrable nonlinear equations, 
Physica 31D (1988), 339-388. 

[12] F.W. Nijhoff, Integrable hierarchies, Lagrangian structures and non-commuting flows, in Top- 
ics in soliton theory and exactly solvable nonlinear equations, (Oberwolfach, 1986), (World Sci. 
Publishing, Singapore, 1987), 150-181. 

[13] F.J. Ernst, New formulation of the axially symmetric gravitational field problem, Phys. Rev. 

167, (1968), 1175-1178. 

[14] F.J. Ernst, New formulation of the axially symmetric gravitational field problem II, Phys. Rev. 

168, (1968), 1415-1417. 

[15] G.A. Alekseev, Soliton configuration of interacting massless fields, Sov. Phys. Dokl. 28 (1983), 
133-135. 

[16] A. Tongas, D. Tsoubelis and P. Xenitidis, A family of integrable nonlinear equations of hyper- 
bolic type, J. Math. Phys. 42 (2001), 5762-5784. 

[17] R. Gamier, Sur des equations differentielles du troisieme ordre dont l'integrale generale est 
uniforme et sur une classe d'equations nouvelles d'ordre superieur, Ann. Ecol. Norm. Sup., 29 
(1912), 1-126. 

[18] J. B. Griffiths, Colliding Plane Waves in General Relativity (Oxford University Press, Oxford, 
1991). 



14 



